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On increasing the coupling strength (A) of a non-Abelian gauge field that induces a generalized 
Rashba spin-orbit interaction, the topology of the Fermi surface of a homogeneous gas of noninter- 
acting fermions of density p ~ k% undergoes a change at a critical value, At ~ &f [Phys. Rev. B 84, 
014512 (2011)]. In this paper we analyze how this phenomenon affects the size and shape of a cloud 
of spin- 1 fermions trapped in a harmonic potential such as those used in cold atom experiments. We 
develop an adiabatic formulation, including the concomitant Pancharatnam-Berry phase effects, for 
the one particle states in the presence of a trapping potential and the gauge field, obtaining approx- 
imate analytical formulae for the energy levels for some high symmetry gauge field configurations 
of interest. An analysis based on the local density approximation reveals that, for a given number 
of particles, the cloud shrinks in a characteristic fashion with increasing A. We explain the physical 
origins of this effect by a study of the stress tensor of the system. For an isotropic harmonic trap, 
the local density approximation predicts a spherical cloud for all gauge field configurations, which 
are anisotropic in general. We show, via a calculation of the cloud shape using exact eigenstates, 
that for certain gauge field configurations there is systematic and observable anisotropy in the cloud 
shape that increases with increasing gauge coupling A. The reasons for this anisotropy are explained 
using the analytical energy levels obtained via the adiabatic approximation. These results should 
be useful in the design of cold atom experiments with fermions in non-Abelian gauge fields. An 
important spin-off of our adiabatic formulation is that it reveals exciting possibilities for the cold- 
atom realization of interesting condensed matter Hamiltonians by using a non-Abelian gauge field 
in conjunction with another potential. In particular, we show that use of a spherical non-Abelian 
gauge field with a harmonic trapping potential produces a spherical geometry quantum hall like 
Hamiltonian in the momentum representation. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85. Lm 



I. INTRODUCTION 

Developments in the area of cold atoms in the past 
decade [IH3] have opened the possibility of obtaining in- 
sights into many outstanding problems of condensed mat- 
ter physics using cold atom quantum emulators. While 
this program has met with early success, serious difficul- 
ties have hampered rapid experimental progress. Notable 
roadblocks are the problem of entropy removal ("cool- 
ing") and the generation of magnetic fields necessary for 
the attainment of interesting states such as the correlated 
quantum hall states. 

There have been many theoretical proposals [1HZ] for 
the generation of both Abelian and non-Abelian gauge 
fields. Interest in this problem has resurged owing to 
recent experimental progress from the NIST group in 
the generation of synthetic gauge fields (both Abelian 
and non-Abelian) using bosonic 87 Rb atoms. [8HT0] A spa- 
tially inhomogeneous Abelian gauge field produces non- 
trivial effects of a magnetic field, while a uniform Abelian 
gauge field is tantamount to a mere gauge transforma- 
tion. On the other hand even a uniform non-Abelian 
gauge field produces interesting physical effects. In the 
context of bosons some aspects of these have been inves- 
tigated and reported. [TMP2] 

A natural question that arises pertains to the effect of 
a uniform non-Abelian gauge field on fermions. Such a 
gauge field induces a generalized Rashba spin-orbit inter- 
action. The hint that such a system contains interesting 



physics came from ref. [13 which demonstrated that a 
high symmetry (see below) non-Abelian gauge field in- 
duces a bound state between two fermions for any attrac- 
tive interaction however small (small negative scattering 
length). This motivated the study [2] of the evolution 
of the many body fermionic ground state with the in- 
creasing strength A of the non-Abelian gauge field. The 
ground state realized for a weak attractive interaction 
in the absence of the gauge field is the BCS (Bardeen- 
Cooper-Schricffer) superfluid state with large overlapping 
Cooper pairs. Ref. [14] demonstrated that on increasing 
the gauge coupling strength A the ground state evolves 
from this just discussed BCS state to a BEC (Bose- 
Einstein condensate) state, i. e., the non-Abelian gauge 
field induces a BCS-BEC crossover even at a fixed weak 
attraction that is unable to produce a bound state be- 
tween two fermions in the absence of the gauge field. 
The most remarkable aspect of the crossover induced 
by the gauge field is that the BEC state obtained for 
large gauge couplings is made up of a new kind of tightly 
bound fermion-pairs called [T4"| "rashbons". Hence, the 
condensate is called as rashbon-BEC (RBEC). Rashbons 
are anisotropic bosons characterized by an anisotropic 
dispersion and spin structure (nematicity), all of which 
are determined solely by the gauge field. The proper- 
ties of RBEC, in particular the transition temperature, 
is determined by those of the corresponding rashbons. 
Properties of rashbons are reported in ref. |15j , where it 
is demonstrated that a non-Abelian gauge field can en- 



2 



hance the exponentially small transition temperature of 
a Fermi gas with a weak attraction to the order of the 
Fermi temperature. Aspects of rashbon superfluidity, ef- 
fects of Zeeman field and effects of imbalance etc. have 
been reported in refs. [TMTU] . 

A non-Abelian gauge field also has quite interesting 
effects on noninteracting fermions. In fact, for a fixed 
density p of the particles, increasing the gauge cou- 
pling strength A induces a change in the topology of 
the Fermi surface at a critical gauge coupling strength 
At- As shown in ref. [H], At is always of the order of 
hp, the Fermi wave vector determined by the density 
(p = fc|/(37r 2 )). In fact, in the presence of a weak at- 
traction, the regime of gauge coupling over which the 
crossover from the BCS to RBEC takes place coincides 
with A > At- In our opinion, one of the first things cold 
atom experiments designed to study fermions in non- 
Abelian gauge fields should probe is this transition in 
the topology of the Fermi surface. The discussion above 
naturally gives rise to many interesting questions. Since 
almost all cold atom experiments are performed in a trap, 
it is useful to know the key signatures that provide mea- 
surable/falsifiable proof of the physics discussed above. 
A specific question is: How the presence of a generic 
non-Abelian gauge field influence the size and shape of 
a cloud of fermions? Answering this question is the key 
motivation of this paper. 

In this paper we investigate noninteracting fermions in 
a parabolic trapping potential characterized by a scale 
cj , in presence of a synthetic non-Abelian gauge field 
(whose strength is characterized by A). Our goal is to 
understand how a synthetic non-Abelian gauge field af- 
fects the size and shape of a cloud of noninteracting 
fermions. To this end: (a) Noting that the most in- 
teresting regime corresponds to luq <C A 2 , we develop 
a Born-Oppenheimer like approximation for the states of 
a trapped particle. This analysis reveals how the inter- 
nal fast degree of freedom induces a Berry connection, 
i. e., a gauge potential on the motion of the particle, 
and most importantly suggests routes to generating in- 
teresting states such as the spherical geometry quantum 
hall states using cold atoms, (b) We study the size and 
shape of the cloud of trapped noninteracting fermions 
under the influence of a non-Abelian gauge field using 
the local density approximation (LDA). (c) Finally, we 
compare the results of the LDA with exact numerical 
calculation. This last study reveals systematic and ob- 
servable deviations from the LDA which are explained 
using the result of our adiabatic theory. This study finds 
that the cloud shrinks (consistent with ref. |17j). Impor- 
tantly we uncover the scaling of the cloud size and the 
evolution of the anisotropy of the density profile with in- 
creasing gauge coupling, characteristic to various gauge 
field configurations of interest. We believe this will be of 
value for experiments on fermions in a non-Abelian gauge 
field. An important byproduct of our adiabatic analysis is 
the possibility of realization of interesting physics in cold 
atomic systems such as the spherical geometry quantum 



hall state by using a synthetic non-Abelian gauge field in 
conjunction with another potential. 

Section [IT] consists of the background including nota- 
tion, estimation of scales and sets up the statement of 
the problem. This is followed by a discussion of the one 
particle states in section |TlT] The effect of the gauge field 



on the cloud size and shape is discussed in section IV 



The final section [V] consists of an itemized summary of 
the paper. Appendix [X] outlines the numerical method 
used to obtain exact one particle states for an extreme 
oblate gauge field configuration. 



II. PREAMBLE AND PROBLEM STATEMENT 

Denoting ^(r) = {^(r)}, a =t,| as creation oper- 
ators of spin-i fermions at position r in three spatial 
dimensions, the Hamiltonian under consideration (see 
refs. [131 Q3] for more details) is 



dr * f (r) ( yl 



Pa 



*(r), 



(1) 



where p is the momentum operator, 1 is the SU(2) iden- 
tity, r M (p — x, y, z) are Pauli matrices, p\ — J2iPi^i e i> 
ei's are the unit vectors in the i-th direction, i — x,y, z. 
The gauge-field configuration (GFC) is described by a 
vector A = J2i^i e i> an d A = |A| is the gauge-coupling 
strength. We work in units where the Planck constant (h) 
and the fermion mass m are set to unity. We are partic- 
ularly interested in high symmetry GFCs called extreme 
oblate (EO) GFC with A = -^(e x + e y ) and spherical 

(S) GFC with A = ^(e K + e y + e 2 ), which nurture in- 
teresting physics. [T3 [ [14 ] 

The eigenstates of the noninteracting Hamiltonian 
(eqn. (FT])) are 



\pa) = \p) <g> \x a (p)} 



(2) 



where \p) is the plane wave state with momentum eigen- 
value p, a — ±1 is the helicity with 



Pa • r\x a (p)) = a\p x \\Xa(p)) ■ 



(3) 



In other words, \Xa(p)) is the eigenstate of p\ ■ r in the 
spin sector with associated helicity a. The energy eigen- 
values associated with the states in eqn. Q are 



£a(p) = y- at\p x \ 



(4) 



We now introduce an isotropic harmonic trapping po- 
tential 



Hn 



w - 



J dr r 2 * t (r)'3>(r), 



(5) 



where ujq is the trapping frequency. Unfortunately, the 
trapping potential spoils the symmetries of the system. 



3 



Not only does it not commute with the usual kinetic en- 
ergy term, it also does not let the helicity a of eqn. ^ 
to be a good quantum number. This renders the diago- 
nalization the full Hamiltonian 



H = Hr+Ht 



(6) 



for a generic ujq and A analytically intractable. 

Fortunately, there are two well separated energy scales 
in the problem in the regime of interest. To see this, the 
typical traps used in experiments are such that the Fermi 
energy (set by the density of particles at the trap center) 
denoted as E F = k F 2 /2 and the trap frequency satisfies 



UJQ 



1(V 



(7) 



We have also noted that the regime of interest of the 
gauge coupling corresponds to A > At where a change 
in the topology of the Fermi surface is engendered by 
the gauge field. In this regime of gauge coupling all the 
occupied states are of positive helicity. Noting that At ~ 
k F for GFCs of interest, we see that the regime of interest 
corresponds to 



LJ 



> io 2 



(8) 



In the next section, we show that the regime indicated in 
eqn. Q allows us to make the adiabatic approximation 
and obtain the energy levels. 



III. ONE PARTICLE STATES 
THEORY 



ADIABATIC 



A. Adiabatic Approximation for a Generic GFC 

The development of the adiabatic approximation to the 
diagonalization of eqn. ^ be gins with the classification 
of the slow and fast degrees of freedom. This is most con- 
veniently done in first quantized notation in momentum 
representation: 



H = — 



if 



2 dp 2 



2 



P\ T 



(9) 



where we have used r — i-g^, the position operator in 
momentum representation. In the regime of interest 
(eqn. Q), A ^> -^/wo, the internal spin degree of free- 
dom is the "fast variable", while the orbital degree of 
freedom is the "slow variable" . This identification of the 
fast and slow degrees of freedom suggests an ansatz for 
the wave function of the particle 



= J dpifj(p) \p) <8> |Xo(p))- 



(10) 



The main physical content of this approximation is that 
the particle motion is such that with the change of p, 



it instantly attains internal spin state corresponding to 
a given helicity a associated with p, hence its internal 
state is \Xa{p))- We shall restrict attention to a = +1 
since this is the relevant case for A > At- 

The goal now is to find the effective Hamiltonian for 
the wave function ip(p). Since the spin Hilbert space has 
nontrivial topology, the associated Pancharatnam-Berry 
phase [213] induces a connection which results in a new 
gauge field and a potential for the slow degree of freedom 
p. Following the line of analysis of ref. [20] , we obtain the 
effective Hamiltonian as 



K 



eff 



up 

2 



d_ 
dp 



e a (p) + V BO (p) (11) 



where 



-i(Xa(p)\ 



dXajp), 

dp 1 



(12) 



is the induced connection (U(l) gauge potential) and 
1 / /dx a (p) l dx a (p) ) 

(13) 



V BO (p) 



2 



dpi dpi 



|Xa(p)XXaCP)|- 



dpi 



dpi 



is a potential where repeated spatial indices are summed. 
The energy levels are now obtained from the eigenvalue 
equation 



U cS ip(p) = eijj(p). 



(14) 



It is quite interesting to note that adiabatic motion in a 
non-Abelian gauge field results in a U(l) (Abelian) gauge 
field for the slow degree of freedom. Clearly, this promises 
to give rise to new possibilities with cold atom systems. 
We shall come back to this point in the sections that 
follow. 

We now discuss application of the formulae developed 
in this section to various GFCs of interest. 



B. Extreme oblate (EO) GFC 



To solve for the levels of eqn. (14) for the EO GFC, 



we choose cylindrical polar coordinates in the momentum 
space (p,<j),p z ) with associated unit vectors e p , and 
e z . Here we consider the positive helicity state (a = +1) 
and 



\X + (P)) 



1 

\/2 V e 



1 

oi<t> 



(15) 



We immediately obtain from eqn. ( 12 1 and eqn. ( 13 1 that 

2 



2p 4 



V BO (p) 



(16) 



The gauge field corresponds to a half flux-quantum mag- 
netic field line running along the z axis of the momentum 



4 



space. The effective Hamiltonian from eqn. (11) reduces 
to 



W 



dp 2 



ld_ 
p dp 



Wo (T 



A) 



v_ 

2 V2 
^o 2 d 2 
2 



A w 2 

-P + 8^ 



(17) 



a P 2 2 



where L z = —i-^r is the z component of the orbital an- 
gular momentum operator. 

We use the separation of variable ansatz ip(p, 4>>Pz) = 
ipp(p)^p^((f>)ip z (z), and note that L z eigenvalue can be 
labeled by an integer to to obtain 



helicity. This approximation is quite valid when the mo- 
mentum of the particle is large (p > A). In the presence 
of the trapping potential, helicity fluctuations are max- 
imum for momenta near the origin. Now note that the 
gauge potential induced by the Berry connection effec- 
tively provides for a nonzero centrifugal barrier for all 
to including m — in momentum space (see eqn. (18 1, 
term proportional to (to + \) 2 )- This effectively "keeps 
the particle away from the origin" in momentum space, 
hence minimizing helicity fluctuations and rendering the 
adiabatic approximation accurate. The same argument 
provides the reason why the approximation works better 
for larger m. In the discussion that follows we shall use 
eqn. (20 1. We note that a similar formula was arrived 
at by different considerations in refs. [23] and [23] that 
appeared during the preparation of this manuscript. 



/ , 2 

W ° "I \ 

— Y u (p>+ 



Up 

2p 2 



pi 

2 



A 

T2 1 



u(p) = iu{p) 
(18) 

^= (a standard substitution), and e is 

the "energy in the (p, </>) degrees of freedom" . This form 
is tailor-made for a semi-classical WKB treatment |21j. 
which gives 



where ippip) 



e(n, to) 



" 4 



A 2 



(19) 



where n is the "radial quantum number" . Remarkably, 
the presence of the potential induced by the connection 
obviates the need for the usual Langer modification [2TI 
122] necessary in semi-classical analysis of this type and 
has the Maslov index of unity at both the turning points. 
Adding in the part of the energy from the z degree of 
freedom we obtain the energy eigenvalues of eqn. (17) 
dependent on three quantum numbers to be 



/ \ ^ 2 / \ k^o 2 

e(n,m,n z ) « — — + (n + n z + l)u + -p- 



(20) 



We now present a comparison of the result of eqn. ( 20 ) 
with the exact numerical solution of eqn. ^ briefly out- 
lined in appendix |Aj Fig. [ija) shows a plot of Ae — 
e(0, to, 0) + ^j- — ujq obtained numerically (solid lines) 
and the same quantity calculated using full numerical 
WKB solution (dashed lines) of eqn. (18). This numeri- 
cal WKB solution very closely matches the approximate 
analytical formula (eqn. (20)) for A > 6 v / wo. We see that 
the full numerical WKB solution has a quite remarkable 
agreement with the exact solution of eqn. ([6| even for as 
small value of gauge coupling as A = 4y / u}o for m = 0, 
and agreement to even lower values of A for higher values 
of m. 

Why does the adiabatic approximation work so well, 
and why does it work better for larger ml This can be 
understood by noting that key ingredient of the adiabatic 
approximation is to "force" the particle to have a fixed 



C. Spherical (S) GFC 

With the reassuring success of the adiabatic approxi- 
mation for the EO GFC, we now turn to spherical GFC. 
The natural coordinate system to treat this GFC via 



eqn. (11 ) is the spherical polar coordinate system in mo- 
mentum space with (p, 6, (j)) as the coordinates and e p , e$ 
and as the basis vectors. The positive helicity eigen- 
states have 



Ix+(p)) 



cos I e-^l 2 
sinfe 1 ^ 2 



It follows from eqn. (12) and eqn. (13) that 



cot 6 
'~2p 



3*, 



V BO (p) 



^0 

4p2 



(21) 



(22) 



Quite interestingly, the gauge potential A corresponds to 
the presence of a monopole of charge Q = \ at the origin 
of the momentum space. The angular motion (i. e., (0, (j))) 
of a particle in a monopole field has been extensively 
studied[2HJ [25] and also used in the context of quantum 
Hall effect in the so called spherical geometry (see section 
3.10 of [27]). We thus obtain, 



H 



cff : 



^0 

2 

2p 2 
Ap 2 



]_d_ 2 d_ 
p 2 dp^ dp 



1 d 

p 2 A 



sm 



89 



Qcote 



d 



sin 9 d<j) 



(23) 



Again using a separation of variable ansatz, ib(p, 9, cf>) = 
tp p (p)Q(9,<f)), and noting that the eigenvalues of the an- 
gular operator in the square bracket in eqn. ( 23 1 is 
1(1 + 1) - Q 2 with t = \Q\,\Q\ + !,••• where I is 
the angular momentum quantum number. For each I 
there are 2£ + 1 states with magnetic quantum numbers 
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FIG. 1. (Color online) One particle levels of a trapped particle in a n EO GFC. Comparison of numerically calculated adiabatic 
WKB energies (dashed lines) obtained from the analysis of eqn. ( 18 1 and the exact result (solid lines in (a) and points in 

I I — \ n 

(b)) obtained using the method outlined in appendi x |A [ Here Ae = e(0,m, 0) + 

indistinguishable with the approximate result of eqn. 



the results are 



! LJ0- -TOr A £ D^/LJO 

(a) Energy for various values of m quantum number as a function 



of A. (b) Energy for A 



/2ujo as a function of the quantum number m. 



rn = —£, — £ + 1, ...,£, and the wave function ft(0, </>) is 
one of the monopole harmonics. [27] We thus obtain 



1 d 2 dtpp 
dp 



P 



2 \p 2 dp* 



(24) 



which is amenable to a semi-classical WKB treatment via 
tppip) — u(p)/p. An approximate analysis gives us 



e{n,£,m) « -— + {n + ~)w + 



(25) 

This result, again, is obtained without the usual Langer 
modification applied to WKB treatment of three dimen- 
sional problems. [22] 

A key inference that can be made from the results 
of this section is that non-Abelian gauge fields used in 
conjunction with a trapping potential can give rise to 
many new possibilities with cold atoms. In particular 
when the gauge coupling is strong compared to the trap- 
ping potential, the motion of the particle is adiabatic 
with respect to the spin degree of freedom. The asso- 
ciated Pancharatnam-Berry phase produces an effective 
gauge field which can be used to realize Hamiltonians 
that are quite interesting. In this regard, it is very inter- 
esting to note that the spherical non-Abelian gauge field 
in conjunction with a harmonic trapping potential just 
discussed produces a "spherical geometry" quantum hall 
like Hamiltonian in the momentum space. Exploitation 



of this aspect of non-Abelian gauge fields should open up 
very interesting possibilities with cold atoms. 

IV. SIZES AND SHAPES OF TRAPPED 
CLOUDS 

In this section we investigate the natural question that 
arises: Is there any discernible effect on the size and 
shape of the trapped clouds due to the influence of the 
non-Abelian gauge field? In the following subsection we 
answer this question using the local density approxima- 
tion and compare it with the exact calculation in the 
subsection that follows. We restrict our attention to zero 
temperature. 



A. Local Density Approximation 

In the LDA [55] , the spatial dependence of the chemical 
potential p is determined by the trapping potential, 



(26) 



where Mo is the chemical potential at the trap center, and 
n(r) is the "local chemical potential" that determines the 
density p{r) of fermions at the point r. It follows from 
eqn. ( p6| that within LDA, fx(r) and hence p(r) will be 
dependent only on r = \r\. The density p is related to 
the chemical potential p, by the equation of state. For 
the GFCs of interest, we can obtain the equation of state 
analytically. We shall focus mainly on the EO GFC for 
which the equation of state is 



() 



,00 = + *fc + + 2 ^) sin " ( J J if ^ °' 

if-¥<„<o. 




X/k\ 



FIG. 2. (Color online) Dependence of the radius of the cloud 
on the gauge coupling A for the EO GFC. Ro is the radius 
of the cloud in the absence of the gauge field and k F is the 
Fermi wave vector related to the trap center density in the 
absence of the gauge field. 



The system is such that the radius of the cloud at 
A = is i?o, with a trap center density given by 
Po = (fc^) 3 /37r 2 and an associated Fermi energy E F = 
(k F ) 2 /2 — ^o 2 Ro- On application of the gauge field, 
the radius of the cloud changes to R. We seek to deter- 
mine the dependence of R on A. This is determined from 
the equation 



47r / dr r p(/j,(r)) 
'o 



24 



RqUJo" 



(27) 



where the right hand side is the number of particles in 
the trap expressed in terms of the radius Rq and the trap 
frequency lo . 

Fig. [2] shows a plot of the dimensionless cloud radius 
R/Ro as a function of the dimensionless gauge coupling 
strength \/k F of the EO gauge field. Quite remarkably, 
the cloud shrinks on application of the gauge field. This 
is consistent with the results of ref. [T7] . Further, we find 
here that there is a critical gauge coupling A c given by 



A,= 



1/6 



k° 

Hp, 



(28) 



such that when A > A c , the radius is given by 

1/5 

for A > A c (29) 



R 
Ro 



5^2 k° F 



This is the regime of gauge coupling that has the most 
interesting physics. This is because for A > A c the local 



Fermi sea every where in the trap will contain only + 
helicity eigen states. In the presence of attractive inter- 
action between fermions in the singlet channel, a rashbon 
BEC is obtained at the center of the trap at low temper- 
atures when A > A c . 

For the S GFC, briefly, the cloud radius satisfies 



R 
Ro 



A 

hp 



R 

Rq 



= I. 



(30) 



Indeed, even for this GFC, the cloud shrinks under the 
influence of the gauge field. As is evident, for X/k F > 1, 

the cloud radius goes as ^ ~ 

The discussion above makes it clear that the shrink- 
ing of the cloud is a characteristic feature obtained on 
application of a generic non-Abelian gauge field. It is 
useful to further investigate the physics behind this re- 
markable and readily observable effect. Let us briefly re- 
cap the factors that determine the cloud size. In a Fermi 
system, the "confining walls" of the container resist the 
Pauli (degeneracy) pressure. In this case equivalent to 
the confining wall is the harmonic trapping potential. It 
is apparent that on application of the gauge field the de- 
generacy pressure reduces, an inference that is consistent 
with the reduction of the cloud size. 

Why should the Pauli pressure of the system reduce 
in the presence of a non-Abelian gauge field? To an- 
swer this question, we consider a homogeneous system 
(no trap) and obtain the stress tensor from the momen- 
tum balance equation. For our noninteracting system, 
the stress tensor operator is (sum over repeated indices 
is implied) 

S ij s ', / dfc*t (fe)fc V CTCT ,(fc)Mv(fc) (31) 
where ^ a (k) is the Fermion creation operator in k space, 
^(fe) = (^-A«rS) (32) 

are the components of the velocity operator (here re- 
peated j is not summed). Note that the velocity opera- 
tor does not commute with the Hamiltonian in eqn. ([!]) . 
A straightforward calculation now shows that the stress 
tensor in our case is 

W=8^ / dfc ^Z kl (v 3 (k)) a n F (s a (k)-p^J 

where n F is the Fermi function, [i is the chemical poten- 
tial and (u J (fe)) a is the expectation value of the velocity 
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We also note that the stress tensor calculated via 



X/k F 

FIG. 3. (Color online) Dependence of the pressure P of a 
homogeneous gas of fermions in an EO gauge field on the 
strength of the gauge coupling A. Po is the pressure of the 
free Fermi gas (no gauge field), and Uf is the Fermi wave 
vector determined by the density. 



operator (eqn. ( 32 ) ) in the helicity eigenstate \ka) given 
in cqn. ©. The pressure P is now obtained as 



For the EO gauge field, 



(v(k)) a = (k 



aX 

71 



)e k + k z e z 



(34) 



(35) 



where we have used cylindrical polar coordinates in k 
space. This formula illustrates a key point. For the 
positive helicity states (a = +1) which are of lower 
energy, the expectation value of velocity vanishes for 
fc = (X/V2, 0, 0). This is the result of the fact that al- 
though the canonical momentum of the low energy states 
is of the order of A, the mechanical momentum which is 
proportional to the velocity is small. It is the mechani- 
cal momentum that contributes to the stress tensor. The 
reason for the fall in pressure is now evident. In the pres- 
ence of the gauge field, all the states with + helicity have 
a smaller velocity (mechanical momentum) . Since in the 
presence of the gauge field, there are always more occu- 
pied states with + helicity, the pressure is expected to 
fall with increasing A for a given density p of particles. 

Fig. [3] illustrates the dependence of the pressure P on 
the gauge coupling A for the EO gauge field calculated 
using eqn. (33) and eqn. (34). The pressure falls as ex- 

pected. In fact, for A > At (At = (3^) k F , the 
gauge coupling at which the Fermi surface undergoes a 
topological transition, k F — [i^p) 1 ^ where p is the 
density), we obtain 



P 
Po 



eqn. (33) is isotropic for this GFC. 



Similar physics is at play also in the spherical gauge 

field where -pr ~ ix) ^ or S au S e couplings larger than 
that required to produce the change in topology of the 
Fermi surface. The more rapid fall of pressure with in- 
creasing gauge coupling owes to the fact that there are 
many more low energy states with vanishing velocity in 
the spherical gauge field. 

These considerations now provide a clear physical pic- 
ture of the shrinking of the cloud with increasing gauge 
coupling. 



B. Cloud Shape 

Within LDA the sole effect of the gauge field is the 
shrinking of the cloud and cloud shape remains spheri- 
cal for this isotropic trap. For the EO gauge field this 
is a somewhat surprising result since one would expect 
the one body wave function to have different behavior 
"in-plane" of the gauge field (x — y plane) and "out of 
plane" (z direction). To investigate this, we calculated 
the density profile using the exact wave functions for this 
problem obtained numerically using the method given in 
appendix |A"| and compared this with the LDA result. 

Fig. [4] shows the density profiles obtained using the 
exact numerical calculation of appendix [X] and that ob- 
tained from LDA for three values of gauge coupling A 
for the EO gauge field. For this gauge field, the symme- 
try of the problem ensures that the density is a function 
only of the in-plane radial coordinate r and the z coor- 
dinate, i. e, p = p(r,z). We show curves corresponding 
to in-plane density profile p(r, 0) and out of plane pro- 
file p(0, z) obtained from the exact numerical calculation. 
The results show a remarkable feature. For A = 8-^/2^0 
we see that in-plane and out of plane density profiles are 
close to each other and agree quite well with the LDA 
result. However, at a larger value of the gauge coupling 
A = 15\/2ujq, 20v / 2wo, we see that the in-plane density 
agrees quite well with the LDA result. However, in both 
of these cases, the out of plane density systematically de- 
viates from the in-plane (and LDA) results. 

This effect can be understood by using the one particle 
spectrum given in eqn. ( |20[ ) obtained using our adiabatic 
theory. Filling up of the levels associated with the quan- 
tum numbers (n,m,n z ) occurs by the ordering of their 
energies. When A is large, the levels of the type (0, m, 0) 

are the lowest lying ones with the spacing between them 

2 

approximately equal to 2m-p-. On the other hand occu- 
pation of a state with n z = 1 or n = 1 requires an energy 
of order ujq higher than the ground state energy. Thus 
the number of states with n = n z = will be determined 
by the maximum value of m attained, m ,nax ~ 2\/^/Uo. 
If N p is the number of particles, we can estimate the 
maximum value of n z of all occupied states. A simple 



for A > Xi 



(36) 



calculation shows that n r ' 



For the small- 
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FIG. 4. (Color online) Density profile of a gas of 910 trapped 
fermions in an EO gauge field. The in-plane density profile is 
indicated by p(r, 0) and out of plane profile by p(0, z). 



est value of A shown in Fig. [4] the value of m max , n" lax 
etc. are all approximately equal and the cloud is approx- 
imately isotropic with reasonable agreement with LDA. 
However, for the higher values of A in Fig. [4j n max sa 6 
but rn max is much larger. The reason why p(r, 0) agrees 
with LDA is now evident. Contributions to p(r, 0) arises 
from many states with different m values and therefore 
the averaging/smoothing of the density is more promi- 
nent. This is not the case for p(0, z) where the density 
contribution is obtained only using n™ ax number of z 
wave functions. Not surprisingly, the out of plane den- 
sity profile p(0, z) shows large deviations from the LDA 
result. Finally, we address the reason of why p(r, 0) falls 
below p(0, z) for larger A. As A increases, the in-plane 
wave functions are "more localized" near the origin. This 
is because the low energy states are predominantly con- 
structed out of plane waves with in-plane momentum of 
order A. This means that real space wave function of the 
particle corresponds to that in which the particle is local- 
ized over a radial distance of order -r . These arguments 
suggest: a) the trap center density (for a given number 
of particles) scales with A, b) the cloud will become more 
and more cigar shaped extending in the z direction as A 



increases. Indeed, our results for the trap center density 
for larger As do agree with this argument, and for very 
high value of A, we do find the tendency for the cloud to 
become cigar shaped. Quite interestingly, the density at 
a given distance from the trap center can be up to 10% 
different along the in-plane and the out of plane direc- 
tions. We believe that this systematic anisotropy should 
be observable and measurable in experiments. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the physics of 
trapped noninteracting fermions in the presence of a non- 
Abelian gauge field that induces a generalized Rashba 
spin-orbit interaction. Specifically: 

1. An adiabatic approximation is developed to ob- 
tain the one particle levels in a trapping potential 
for a generic gauge field. The effective adiabatic 
Hamiltonian for the particle includes a gauge po- 
tential arising due to the induced connection from 
the internal spin degree of freedom. For the spher- 
ical gauge field this gauge potential is shown to be 
equivalent to the field of a monopole at the origin 
in the momentum space. Approximate analytic for- 
mulae for the levels are obtained for high symmetry 
gauge field configurations of interest. For the ex- 
treme oblate GFC these results are compared with 
the exact results obtained numerically and excel- 
lent agreement is demonstrated. 

2. The effect of the gauge field on the size of the 
trapped cloud of fermions is investigated within 
a local density approximation. It is shown, quite 
generically, that the cloud shrinks with the increas- 
ing strength of the gauge field. Formulae for the 
dependence of the cloud size on the strength of the 
gauge field are obtained. The physical origins of 
this phenomenon are elucidated by our analysis of 
the stress tensor. 

3. The density profile obtained using the local den- 
sity approximation is compared with that obtained 
from a exact numerical calculation. It is shown 
that for the EO GFC, the density distribution is 
anisotropic in a characteristic way and should be 
observable in experiments. 

We believe that the results obtained in the paper will 
be useful to design experiments involving fermions in 
non-Abelian gauge field. In addition, our results also 
open doors to new possibilities with cold atoms that 
arises from using the non-Abelian gauge field in conjunc- 
tion with a trapping potential or other types of poten- 
tials. In particular, we have shown that a spherical non- 
Abelian gauge field and a harmonic trapping potential 
can be used to simulate the spherical geometry quantum 
hall system. Interestingly, the quantum hall system in 
this case is realized in the momentum space. 
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Appendix A: Numerical Calculation of One Particle 
States for the EO GFC 



In this appendix we outline the method used to nu- 
merically obtain the exact one particle states for trapped 
fermions in an EO gauge held. 

As is evident we need to focus only on the in-plane 
degrees of freedom which are most conveniently de- 
scribed using polar coordinates (r, 4>) . The Hamiltonain 
(eqn. Q) is rewritten in a more convenient form 



H = 



w 



2 



V2 



P + T 



P T + ) 



(Al) 



H 



where p^ 1 = p x ±ip y and r = s(t x ±ir y ). In the follow- 
ing, energy will be measured in the units of loq, distance 
in units of 1 /y/u>o, momentum and gauge coupling in the 
units of Jwo, rendering all quantities dimensionless. 
The eigenstates of Ho are such that 

U \m,n,a) = (2n + |m| + l)\m,n,a) (A2) 

where m is the angular momentum quantum number and 
n is the radial quantum number, and a is the spin state 



which is j" or 4-, i.e., 

|m, n, a) — \m, n) ® \a) 

with 



(A3) 



(r,(j)\m,n) = 



n{n + m)\ 



> m L™(r 2 )e lm ^ (A4) 



where L™ are associated Laguerre polynomials. 
Since J-, = L, 



jr z commutes with the Hamitlonian 
(eqn. (Al)), we construct energy eigenstates which are 
also eigenstates of J z . Noting that state |m, n, f) and 
|m + l,n, J,) are eigenstates of J z with the eigen value 
m+h, we can write the eigenstates of Hamiltonian ( Al ) 
as 



1$) = 



En 6 « 



„„|m,n,t) 
|m + 1, n, J.) 



(A5) 



where a mn and b mn are coefficients to be determined. 
This leads to a matrix eigenvalue problem which is com- 
pletely defined with the specification of the following ma- 
trix elements: 



(m + l,n,i \H\\m,n',t) = (Vm + n + 1 S n<n ' +Vn+ 15„,„'_i) 

(m,n,t \H\\m + l,n',l) = (Vm + ri + 1 8 n , n > + V?i' + 1 ^ n ',n-i) 



(A6) 



In actual calculations the sum in eqn. (A5| has to be 
truncated at an upper value iV of n. We have chosen N in 



such a way that the energy eigenvalues are converged to 
the desired double precision accuracy. When A = 20V%, 
we used N = 180 to obtain this accuracy. 
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